New universality class for the fragmentation of plastic materials 
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We present an experimental and theoretical study of the fragmentation of polymeric materials 
by impacting polypropylene particles of spherical shape against a hard wall. Experiments reveal 
a power law mass distribution of fragments with an exponent close to 1.2, which is significantly 
different from the known exponents of three-dimensional bulk materials. A 3D discrete element 
model is introduced which reproduces both the large permanent deformation of the polymer during 
impact, and the novel value of the mass distribution exponent. We demonstrate that the dominance 
of shear in the crack formation and the plastic response of the material are the key features which 
give rise to the emergence of the novel universality class of fragmentation phenomena. 

PACS numbers: 62.20. Mk; 46.50.+a; 64.60.-i 



Fragmentation phenomena are ubiquitous in nature 
and play a crucial role in numerous industrial processes 
related to mining and ore processing [l|. A large vari- 
ety of measurements starting from the breakup of heavy 
nuclei through the usage of explosives in mining or frag- 
menting asteroids revealed the existence of a striking 
universality in fragmentation phenomena [H4lol|: frag- 
ment mass distributions exhibit a power law decay, in- 
dependent on the type of energy input (impact, explo- 
sion, ...), the relevant length scales or the dominating 
microscopic interactions involved. Detailed laboratory 
experiments on the breakup of disordered solids have re- 
vealed that mainly the effective dimensionality of the sys- 
tem determines the value of the exponent, according to 
which universality classes of fragmentation phenomena 
can be distinguished. Several possible mechanisms have 
been put forward to understand the emergence of the 
universal power law behavior. For rapid break-up of het- 
erogeneous bulk solids with a high degree of brittleness, 
the self-similar branching-merging scenario of propagat- 
ing unstable cracks governed by tensile stresses can ex- 
plain the main features of the fragment mass distribution 
[5l. [ll| - [l3| . while for shell systems an additional sequential 
binary breakup mechanism has to be taken into account 
0, El • It is an important question of broad scientific and 
technological interest how plasticity, and the emergence 
of complicated stress states like shear affect the fragmen- 
tation process. The fundamental questions of how robust 
the universality classes are with respect to mechanical 
properties and whether there exist further universality 
classes of fragmentation of solids, still remain open. 

In the present Letter we investigate the fragmentation 
process of plastic materials by impacting spherical parti- 
cles made of polypropylene (PP) against a hard wall. Our 
experiments show that the mass distribution of plastic 
fragments exhibits a power law behavior with an expo- 
nent close to 1.2, which is substantially different from the 




FIG. 1: (Color online) Final states of impact at low impact 
velocities in the experiment (a) and in the simulation (6). In 
the contact area with the hard wall large permanent deforma- 
tions occur due to compression, while above it vertical cracks 
are formed due to tensile stresses. The simulations are in very 
good agreement with the measurements. 

one of bulk brittle materials in three-dimensions. In or- 
der to understand the physical origin of the low exponent, 
a three-dimensional discrete element model is developed 
where the sample is discretized in terms of spherical par- 
ticles connected by elastic beams. To capture the frac- 
ture mechanisms of plastic materials in the model, broken 
particle contacts are able to reconnect when compressed 
against each other leading eventually to the healing of 
cracks and to a plastic energy dissipation. By adding 
two novel ingredients into a DEM model, namely healing 
of broken bonds and breaking under compression, we are 
able to simulate the plastic behavior of PP and reproduce 
in fact the observed new power law distribution. 

In the experiments we used isotactic polypropylene, 
which is a thermoplastic polymer with an intermediate 
level of crystallinity in its molecular structure. The most 
important parameters of PP are: Young modulus 1300 
MP a (room temperature), glass transition temperature 
-10 °C, melting point 160 °C, and density 0.9 g/cm 3 . 
Besides tacticity, the mechanical response and fracture 
characteristics of PP are also strongly affected both by 
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the temperature and by the rate of loading: increasing 
strain rate gives rise to a more brittle response while the 
raising temperature enhances ductility [141 ]. To achieve 
fragmentation a single particle comminution device was 
used which accelerates particles one-by-one by centrifu- 
gal force in a rotor up to the desired velocity. The rotor 
ensures that the particles hit the hard wall at a rectan- 
gular angle in an evacuated environment eliminating the 
disturbing effect of inclined impact and of turbulent air 
flow. Spherical PP particles of diameter d = 4 mm were 
fragmented at different impact velocities vq in the range 
30 m/s-180 m/s. Figure HJa) shows that at low enough 
velocities, the collision does not result in a breakup, in- 
stead the particles undergo a large plastic deformation at 
the impact site. Above the completely flattened contact 
zone of permanent deformation meridional cracks form 
due to tensile stresses, however, the body does not fall 
apart, the plastic zone retains its integrity. 

Fragmentation occurs when the impact velocity i>o ex- 
ceeds a material dependent critical value v c , which is 
about 60m/s for our PP particles. At each impact ve- 
locity 400 particles were fragmented accumulating the 
fragments in the grinding chamber of the machine. In 
the data analysis, 99 — 99.5% of the total mass of the 
samples was recovered. In order to evaluate the mass dis- 
tribution of the fragments, we scanned the pieces with an 
open scanner obtaining digital images where fragments 
appear as white spots on the black background @, @|- 
This way the identification of fragments is reduced to 
cluster searching of white pixels. The very fine powder 
of extension smaller then the pixel size of the scanner 
was left out of the data analysis. The two-dimensional 
projected area w of fragments is determined as the num- 
ber of pixels of the clusters, from which the mass m of 
fragments can be estimated asm~ ui 3 / 2 since the three- 
dimensional fragment shape is close to isotropic. We did 
check the shape isotropy by calculating the square root 
of the ratio of the larger I\ and smaller I2 eigenvalues 
of the tensor of inertia of the 2D projections. The inset 
of Fig. [5] shows that the value of \fl\jl% is close to 1.9 
for almost all fragment masses indicating a high level of 
isotropy. The mass distribution F(m) of fragments ob- 
tained at three different impact velocities is presented in 
Fig.[2j It can be observed in the figure that at the highest 
impact velocity vq = 75 m/s, where the state of complete 
breakup is reached, a power law functional form emerges 
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over more than 3 orders of magnitude in the regime of 
small fragments. At lower impact velocities the power 
law regime of the distribution is followed by a hump for 
the largest fragments which gradually disappears and the 
cutoff becomes exponential as vq increases. The most 
astonishing feature of the experimental results is that 
the value of the exponent t p i = 1.2 ± 0.06 of the power 
law regime is significantly lower than the values ~ 
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FIG. 2: (Color online) Mass distribution of fragments ob- 
tained from experiments at three different impact velocities, 
mo denotes the average mass of PP spheres. At the highest 
impact velocity vo = 75 m/s, F(m) shows a power law be- 
havior over 3 orders of magnitude followed by an exponential 
cutoff for the very large pieces. Simulation results obtained 
with the parameters Qth = 1 and th = (bending&healing) 
are in a very good agreement with the experimental findings. 
Inset: the shape parameter \JhJl2 of fragments as a function 
of mass for the experiments of the main figure. 

1.8 — 2.1 typically found in the fragmentation of three- 
dimensional bulk objects consisting of disordered brittle 
materials EH, Hif . The anomalously low value of t p i 
is the consequence of the breakup mechanism of plastic 
materials which has not been considered by the usual 
theoretical approaches 0, [l2j ■ 

In order to reveal the underlying physical mechanisms 
of the fragmentation of plastic materials, we used a Dis- 
crete Element Model (DEM) to simulate the fragmenta- 
tion of polymeric particles of spherical shape when they 
impact a hard wall. The spherical sample is represented 
as a random packing of spheres with a narrow size dis- 
tribution. The particles are connected by beams along 
the edges of a Dclaunay triangulation of the initial par- 
ticle positions such that in the initial state neighboring 
elements are in contact. In 3D the total deformation of 
a beam is calculated by the superposition of elongation, 
torsion, as well as bending and shearing in two different 
planes [15| . To capture one crucial aspect of the deforma- 
tion behavior and the fracture of plastic materials, our 
DEM model has two novel components, i.e. the form of 
the breaking criterion of the beam elements and the re- 
activation of broken contacts under pressure. The beams 
break when they get overstressed according to a physical 
breaking rule 
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where e denotes the longitudinal strain, 6i and Qj arc the 
generalized bending angles at the two beam ends, and Eth 
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FIG. 3: (Color online) (a) The average mass of the largest 
{M m ax} and second largest (M^^x) fragments as a function of 
the impact velocity «o normalized by the critical velocity v c of 
fragmentation for two parameter sets. (6) Fragment mass dis- 
tributions obtained for the limiting cases of tension, bending, 
and bending&healing dominated breakups giving the corre- 
sponding exponents Ttr, T s h, and T s h-pi, respectively. 



and Q t h are breaking thresholds which have fixed values 
for all the beams. The two terms of Eq. ([2j characterize 
the contributions of the stretching and bending failure 
modes of a beam. Simulations showed that the local 
shear of the particle contacts provides the main contri- 
bution to the bending angles Oi, Qj, so that bending 
dominated beam breaking in Eq. ([2]) characterizes crack 
formation due to shear. Varying the values of the break- 
ing thresholds Eth and Qth, the relative importance of 
stretching and bending can be controlled: increasing the 
value of a breaking parameter, the effect of the corre- 
sponding failure mode diminishes. An important feature 
of the breaking criterion is that in Eq. ([2]) the deformation 
e is not restricted to positive values contrary to former 
simulation studies on brittle fracture and fragmentation 
@, B G3 ■ Since the first term of Eq. (|2|) becomes nega- 
tive when the beam is compressed, failure is dominated 
by the bending/shear mode in such a way that increasing 
compression increases the shear resistance of the beam. 
In order to represent the plastic behavior of the material, 
we assume that the beams have a linear elastic behavior 
up to fracture, but, whenever two particles are pressed 
against one another for a time longer than th, a new, 
undeformed beam is inserted between them. This way 
during the impact process, the particle contacts may un- 
dergo a sequence of breaking-healing events which leads 
to plastic energy dissipation and to the appearance of 
permanent deformation. Varying the healing time th the 
mechanical response of the model material can be con- 
trolled: th — corresponds to the case of perfect shear 
plasticity, while th — >• oo implies no healing at all, i.e. 
brittle behavior. 

In order to investigate the effect of local failure modes 
of beams on the fragmentation process, we carried out 
computer simulations by setting the stretching thresh- 
old to a fixed value e t h = 0.02 and varying the bending 
threshold within a broad range 1.0 < Qth < 200. In this 
way Qth = 200 and Q t h = 1 imply total tension and 
bending dominance, respectively, while intermediate Qth 
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FIG. 4: (Color online) Probability distribution of the angle $ 
between the velocity vector of fragments and the direction of 
impact for tension Q t h = 200 and shear Q t h = 1.0 dominated 
breakup. Insets: final states of fragmentation for tension (a) 
and shear (b) dominance. The arrows indicate the direction 
of the velocity of a few fragments. 



values interpolate between the two limits. Simulations 
were stopped when the system reached a relaxed state, 
i.e. no beam breaking occurred during 1000 consecutive 
time steps. Figure presents the final state of an im- 
pact simulation obtained with a sample of N = 24000 
particles at low impact velocity. It can be observed that 
the model is able to reproduce both the deformation and 
the crack structure of PP, with parameter values where 
the beam breaking is dominated by bending Q t h — 1-0, 
furthermore, compressed contacts easily heal th = 0. The 
large permanent deformation of the sphere in Fig. HJ6) 
arises due to breaking-healing sequences of particle con- 
tacts in the compressed impact zone. Above this zone 
tensile stresses arise resulting in opening cracks along the 
impact direction in agreement with the experiments. Fig- 
ure [3ja) presents the average mass of the largest (M max ) 
and second largest (M^f x ) fragments as a function of 
the impact velocity Vq for the tension Q t h = 200 dom- 
inated brittle breakup th = oo (tension), and for the 
bending dominated fragmentation Q t h = 1-0 with per- 
fect shear plasticity th = (bending&healing). For low 
impact velocities the sample gets only damaged, hence, 
the largest fragment comprises nearly the entire mass Mq 
of the body (M max ) ~ Mq, while the second largest one is 
orders of mag nitude smaller (M* a d x ) < (M max ). Frag- 
mentation occurs when the largest and second largest 
pieces become comparable, i.e. at the critical velocity v c 
where the (M^f x ) curve has a maximum coinciding with 
the inflexion point of (M max ). Figure [3fa) shows that 
for the breakup of brittle materials dominated by tensile 
stresses, the damage-fragmentation transition is sharp in 
agreement with experiments P.[To|. however, when shear 
breaking dominates a broad critical regime emerges. The 
reason is that at v c the largest fragment does not break 
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up into large pieces as for the tension case, but instead 
gradually erodes with a cleavage mechanism giving rise 
to a slowly decaying residue. 

Representative examples for the mass distribution of 
fragments are presented in Fig. [3^6) for three parameter 
sets. For the mass distribution of the tension dominated 
breakup of heterogeneous brittle materials obtained at 
the parameter values ®th = 200, th = oo {tension) a 
power law behavior is evidenced with the usual expo- 
nent Tj, r = 1.9 ± 0.1 [l5|. Simulations showed that de- 
creasing the healing time th in the tension limit 
practically does not affect the fragmentation process be- 
cause fragments are only generated by opening cracks 
which do not let healing play any role. It is impor- 
tant to emphasize that the breakup process substan- 
tially changes when shear dominates the crack formation 
©t/i = 1 {bending). Even in the case of perfectly brittle 
beam breaking th = oo the low shear resistance leads to 
a significantly lower exponent than in the tensile limit 
r s h = 1.0 ± 0.05 (see Fig. Increasing the strength 

of plasticity th — > when shear dominates the cracking 
<d t h = 1, the exponent of the fragment mass distribution 
slightly increases: due to the healing of cracks fragments 
can merge which decreases the relative frequency of large 
pieces leading to a faster decay of the distribution. The 
highest value of the exponent T s h- P i — 1.25 ± 0.06 is ob- 
tained in the plastic limit th — together with Q t h = 1-0 
{bending& healing). It is important to emphasize that 
varying solely the velocity of impact in Fig. [2 for the 
mass distribution of fragments an excellent agreement is 
obtained between the shear-plastic simulations and the 
experiments. The results show that the shear domi- 
nated cracking together with the healing mechanism of 
compressed crack surfaces are responsible for the unique 
fragmentation of plastic materials. During the impact 
process of the experiments a considerable fraction of the 
kinetic energy of the particle gets transformed into heat 
which enhances the effect of healing and plasticity in con- 
sistence with our theoretical results. Our extensive sim- 
ulations indicate that the exponent r s h- p i is universal, 
i.e. it does not depend on the impact velocity or on ma- 
terials' microstructure characterizing a novel universality 
class of fragmentation phenomena. 

Due to low shear resistance, we find that the fragmen- 
tation of plastic shows similarities to the breakup of liq- 
uid droplets colliding with a wall. The spatial distribu- 
tion of fragments and the crack structure in the final state 
of brittle and plastic fragmentation are compared in the 
insets of Fig. |4] Although meridional and segmentation 
cracks are clearly observed in the brittle case in agree- 
ment with experiments, a lateral spreading of fragments 
is obtained for plastic with shear dominated breaking. 
For the quantitative characterization of the spatial spread 
of fragments, the distribution of the angle 4> of the veloc- 
ity vector of fragments with the impact direction (verti- 
cal) is presented in Fig. @] It can be observed that most 



of the fragments bounce back from the hard wall when 
tension dominates, however, in the shear-plastic case the 
fragments escape laterally producing a "splash" of the 
entire body similar to liquid droplets 13, [l8( . 

In conclusion, our experimental and theoretical study 
revealed that the breakup of plastic materials falls into 
a novel universality class of fragmentation phenomena 
characterized by the new mass distribution exponent. 
Based on discrete element simulations we showed that 
the plastic behavior of the material together with the 
dominance of shear in crack formation are responsible 
for the substantial difference from brittle fragmentation. 
The low shear resistance of the material gives rise to a 
splashing similar to the breakup of droplets of highly vis- 
cous liquids. Beyond the industrial importance of the 
fragmentation of polymeric materials, our results might 
also be applied to obtain a deeper understanding of the 
fragmentation of highly viscous magma during pyroclas- 
tic activity at volcanic eruption |l9l I20I]. For theoretical 
investigations our results demonstrate that the breakup 
of solids cannot be understood as a generic stochastic 
process since the precise mechanism of crack initiation 
and growth, i.e. the dominance of tensile or shear stresses 
govern the exponent of fragmentation. Still remains the 
challenge to construct a theoretical approach which ex- 
plains the emergence of universality based on crack dy- 
namics. 
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